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INVOLUTE STRIP DEVELOPMENT METHOD FOR FABRICATION OF ROTATIONALLY 
SYMMETRICAL DOUBLE CURVED SURFACES OF REVOLUTION 
by Rene E. Chambellan and Eugene J. Pleban 
Lewis Research Center 

SUMMARY 

Many methods are available for the approximate development of double curved sur- 
faces of revolution onto a plane. The method presented in this report differs from 
others in that a set of constant-width curved strips (ribbons) or tubes can be used to re- 
produce any surface of revolution. This feature is of particular interest in the design of 
heat exchangers, where a miiform distribution of heat-transfer tubes may be required 
over a surface of revolution. Two important examples of this type are tube-type rocket 
nozzles and fuel-cooled nozzle plugs for jet engines. 

In application to a tube-type heat exchanger, the tubes are preformed into a curve 
matching the elemental curved strip from the shell development. The tube-to- 
elemental- strip joint is easily made and inspected. After joining, the tube-strip assem- 
bly can be formed over a mandrel and joined at the strip edges to form the shell. The 
tubes can be either internal or external to the shell. 

The surface development is performed in two steps. First, the surface of revolu- 
tion is divided into a series of coaxial intersecting right- circular cones, cylinders, and 
annulii, as appropriate. These single curved surfaces are then developed onto a single 
plane where each of these surface elements must have a particular orientation with re- 
spect to its contiguous element. Second, these surface elements are then further sub- 
divided into curved strips of constant width. Winding these strips over a mandrel, rep- 
resenting the original surface of revolution, produces the surface of revolution. The 
development can be accomplished either graphically or by means of a digital computer 
program; both methods are presented in this report. 

The development procedures presented in this report are limited to rotationally 
symmetrical sxorfaces of revolution. However, extension of the principle can be made to 
include more complicated surfaces. It is not possible to completely cover, with 
constant- width strips with no overlapping, closed surfaces of revolution such as spheres 
and streamlined bodies. The curved strips have a finite width and therefore require a 
finite pole circle whose circumference is equal to the strip width times the number of 
strips. An exception to this is the development of a toroid, since it is possible in this 
case to completely cover the surface. 



INTRODUCTION 


This report presents a method for developing a double curved surface of revolution 
into flat ribbons. The development method is presented both as a graphical procedure 
and as a digital computer program. The development is accomplished by developing a 
series of cones, cylinders, and plane annuUi, which are used to closely approximate the 
original surface, and by then subdividing the developed surfaces into curved ribbons. 

Many methods are available for producing double curved thin shells of revolution 
(refs. 1 and 2). Generally, the material for the formation of a double curved shell of 
revolution is initially in the form of flat sheet stock, conical or cylindrical shells, small 
preformed shell segments, or thick shells As appropriate, these material forms can 
be spun, stretch-formed, bent, machined, and welded to produce the desired shape. 

In the course of developing a construction method for a fuel-cooled jet engine nozzle 
plug, where the cooling tubes are to be located and uniformly spaced on the inside sur- 
face of the plug shell, it became apparent that direct access to the tube-shell joint would 
be necessary to ensure an adequate connection. A method, described in this report, has 
been evolved whereby a rotationally symmetric shell can be developed approximately 
into identical planar ribbon strip-elements. In the case of a nozzle plug, these strips 
are joined with the cooling tubes, and then these subassemblies are formed to produce 
the plug shell. The identical procedure can be used to produce subassemblies for the 
construction of tube-t3^e rocket nozzles. 

The development method presented herein is general and not restricted to heat- 
exchanger -type structures. Surfaces where the axis of symmetry pierces the surface of 
revolution reqviire a pole circle of finite size. The holes in the surface of revolution 
represented by the pole circles can be closed, in practice, by suitably formed shallow 
double curved surfaces. A toroidal surface does not have this limitation. 

SYMBOLS 


A 

a, It., . . 

f 

h 

N 

n 

r 


area 

geometrical construction points 
operator, dependent on stated conditions 

projected length of element of surface of revolution on axis of symmetry 
number of ribbons 
ordinal nvunber 

radius of surface of revolution about axis 
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II 

5 curvilinear length of composite involute 
u, V, w, X, y orthogonal coordinate axes 

y angle between radius vector and axis normal to evolute where involute 

starts 

A increment 

rj angle between tangent to an involute and axis normal to evolute where in- 

volute starts 

6 angle about u-axis measured in plane parallel to v-w plane 

v number of elemental composite involute strips 

p slant height of conical element or radius vector in plane of development 

(p one-half vertex angle of right- circular cone 

increment in the angle 9 

CO width of elemental composite involute strip 

Subscripts: 


i index number, 1 or 2 

j index number, 1 or 2, identifying extremities of element of surface of rev- 

olution 

m integer, identif 5 dng intermediate radii of cone n 

n integer identifying particular element of subdivided surface of revolution 

s identifies locations between j = 1 and j = 2, 1 ^ s < 2 

0 indicates involute generating circle 


SUBDIVISION OF SURFACES OF REVOLUTION INTO INVOLUTE STRIPS 
Subdivision of Plane Surfaces 

The development of right- circular truncated cones and cylinders onto a plane results 
in geometric figures which are segments of annulii or rectangles, respectively. Plane 
annulii and segments of annulii can be subdivided by families of involute curves whose 
evolutes (circles herein) are concentric and equal to or smaller than the inner circle of 
the annulus considered. K the spacing of the involutes in the family is made uniform, a 
plane annulus can be subdivided into a number of identical curved strips whose sides are 
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Figure 1. - Plane circular annulus subdivided by involute strips. 
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Figure 2. - Rectangle subdivided by involute strips, with 
center of evolute circle at infinity. 


parallel (see ref. 1), as shown in figure 1. Recognizing that the cylinder is the limiting 
case of a cone (with zero apex angle), then the development will be a rectangle and the 
family of involute curves will degenerate into straight parallel lines, as shown in fig- 
ure 2. These straight lines are developed helices and are the limiting case for the in- 
volute, since the developed helix can be generated by an involute generating circle lo- 
cated at infinity. 


Subdivision of Double Curved Surfaces of Revolution 

Double curved surfaces of revolution cannot be exactly developed onto a plane. 

These surfaces can be approximated by a series of colinear right- circular cones, right- 
circular cylinders, and plane annulii, as appropriate. Since these latter surfaces are 
single curved, they can be developed. The larger the number of surface-approximating 
elements, the greater the accuracy of the development. 

A typical example of the development of a double curved surface of revolution into 
curved ribbons is shown in figure 3. The top of the figure shows a 15°-15° isometric 
view of the surface, and the lower portion of the figure shows the developed curved 
ribbon. 
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Figure 3. - Open-ended streamlined surface of revolution (jet engine nozzle plug). 


Development of Typical Surfaces of Revolution 

In the development of surfaces of revolution such as partial spheres, annulii, rocket 
thrust chambers, and open-ended streamlined bodies, the smallest pole circle that can 
be covered by the strips, with no overlapping, will have a circumference equal to the 
product of the strip widths and the number of strips. Examples of these representative 
bodies of revolution are shown in figures 3, 4, 5, and 6 which are an open-ended 
streamlined surface, an open-ended sphere, a toroid, and a rocket thrust chamber, re- 
spectively. These figures show an outline of the body, the trace of one involute ribbon 
strip centerline and the developed involute strip. These figures were drawn from data 
generated by a computer program described in appendix B based on equations derived in 
appendix A. 
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Centerline of composite 
involute ribbon ^ 



15° - 15° Isometric projection 



Figure 4. - Open-ended spherical surface. 



(entire ribbon is series of these segments). 


Figure 5. - Toroid. 
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Centerline of composite 




Developed composite involute ribbon 


Figure 6. - Rocket thrust chamber assembly. 

Accuracy of Method 

The accuracy of the surface of revolution constructed from the developed involute 
ribbon strips is dependent on the strip width and the number of single curved surface- 
approximating elements such as cones, cylinders, and plane annulii. The greater the 
number of strips and elemental surfaces of revolution, the more accurate the develop- 
ment will be. Thus, in practice, the original double curved surface of revolution can 
be approximated as closely as desired. 

The precision of the development procedure is greatly enhanced when a numerical 
computation method is used in place of a graphical method. Graphical methods are use- 
ful for checking results and for preliminary design; however, a final design should be 
done by the numerical computation method. The precision of the graphical method where 
very small cone apex angles are encountered suffers because very large radii of curva- 
ture are required in the plane development and the involute generating circle is far re- 
moved, in many cases far off the drawing board. In either method of development, 
errors in the development process are cumulative; however, in the numerical method, 
these errors can be reduced by increasing the precision of the numbers used. In the 
graphical method, the only way to improve precision is to enlarge the scale of the draw- 
ing, which quickly becomes impractical. 


Method of Development 

The single curved surface of revolution approximating elements can be established 
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so as to obtain an equivalent surface where either the external or the internal envelope 
is preserved. In either case, the axis of symmetry is subdivided into lengths which are 
not necessarily equal. In the former case (external envelope), the radii to the actual 
surface are obtained at the selected points along the axis. In the latter case (internal 
envelope), the radii to the intersections of contiguous approximating elements, which are 
tangent to the actual surface, are obtained. The relative magnitudes and positions of 
consecutive radii establish whether the approximating element is a cone, a cylinder, or 
an annulus. If a mathematical formula is available, describing the shape of the surface, 
the former approximation is less cumbersome, assuming that the determination of radii 
is a computational instead of a graphical procedure. Neither approach is superior for 
graphical measurements. 

The development of surfaces of revolution can be more lucidly explained as a graph- 
ical procedure. Expressing the procedure in terms of mathematical symbolism is then 
more easily accomplished by the use of geometric relations. 

An open-ended streamlined surface of revolution has been selected as an example 
and developed by a graphical method. Figure 7 (a) shows a longitudinal cross section of 
the surface of revolution. The actual surface is indicated by a dot-dash line. Truncated 
cones, shown as solid lines in the cross section, are used to approximate the original 
surface. These cones are the sirnface-approximating elements and are shown as seg- 
ments of plane annulii in the plane of development (fig. 7(b)). The relative positions of 
the developed cones are dependent on the conditions which require that the segments of 
the composite involute curve be continuous and tangent at the juncture of adjacent devel- 
oped cones. The developed cones must also be tangent at the junction point. 

The reason for using involute curves for the subdivision of the developed surface- 
approximating elements is that a family of involutes are all parallel (ref. 3). Thus, a 
plane surface can be subdivided into any arbitrary number of similar contiguous areas. 
The use of many involutes in the family produces subdivisions which are ribbon- like in 
configuration. 

Since the surface -approximating elements (cones in our example) are single curved, 
their development is exact and, therefore, any subdivision of the development will suffer 
no distortion when the cones are transformed to surfaces of revolution. This conclusion 
leads to the concept that the conical siirface can also be constructed of strips defined by 
a family of involutes in the development plane. 

A complication arises when the involute strips are to be made continuous (in the de- 
velopment plane) in going from one surface-approximating element (cones) to the next. 
Since the developed contiguous surface elements can be tangent at a single point only on 
each boundary, the edges of the involute strip in proceeding from one element to the next 
would suffer a discontinuity. However, a smooth and continuous, single composite- 
involute curve can be defined which passes through each juncture point of the developed 
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(a) Half cross secfion of open-ended streamlined surface of rewlution showing surface-approximating 
cones where pjj^ PlO^ P 52 - 



(bl Graphical development of double curved surface of revolution shown in figure 7(a) (see fig. 10 for method of generating invoiute). 

Figure 7. - Graphical development of open-ended streamlined surface of revolution. (Circles designated by pno = 1,2,3,...) 
are involute generating circles.) 


surface -approximating elements. This composite curve is then used as the centerline of 
an involute ribbon strip and to locate all the points of tangency of the developed surface - 
approximating elements. 

In this example, the development is started with cone 1, which is shown as a segment 
of a circular annulus in the plane of development. The radius of the involute generating 
circle is determined from an assumed number and width of strips, consistent with the re- 
quirement that the product of number of strips times strip width cannot exceed the small- 
est circumference of the surface of revolution. In appendix A it is shown that, in general, 
PnO = 2 ^ ^in y ' therefore, the radius of the first involute general circle is 


n 


^10 


10 


sm 


<Pl 


where r jq = v(x>/2n (from eq. (A7)) and is limited by the requirement that r^j^ ^ r^^g and 
^11 = ^ 10 - 
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The center of each involute generating circle is coincident with the arc center of its 
corresponding developed cone. The first involute curve is drawn, in the plane of develop- 
ment, passing through the inner cone arc on to the exterior arc of the cone where it is 
terminated. This termination point is where the developed cone 2 is then made tangent 
to cone 1. This means that the center of the arc for developed cone 2 is colinear with a 
line drawn through the point of cone tangency and the center of cone 1. In figure 7 (b) 
this line is indicated as ^ 0 . 

At the arc center for cone 2, a new involute generating circle is drawn whose rad ius 
is P 2 Q. Graphically, this radius is obtained b y ex tending the involute generating line ^ a 
to o . Point c is where a line perpendicular t o h a , extended, passes through the arc 
center of developed cone 2. The line segment e o is then the radius pgQ, obtained from 
equation (A9). 

The point It- on the li ne ^ c traces the involute of circle e on the developed cone 2, as 
the extended tangent line It c rolls around circle e. This involute segment is generated 
until it intersects the outer arc of developed cone 2. This process is repeated for suc- 
cessive cones. The arc of developed cone 4 has its center on the other side of the com- 
posite involute curve for cones 1, 2, and 3; which is due to the fact that the slope of 
cone 4, in figure 7(a), is opposite to that for cones 1, 2, and 3. In general, this alter- 
nation of the developed cone centers will occur whenever the composite involute curve 
crosses its tangent. Using this composite involute curve as a centerline, a ribbon of 
width o) can now be drawn which is a developed surface element of the double curved 
surface of revolution. 

Two deviations from the original surface of revolution occur in the surface as con- 
structed from the developed involute ribbon strips. These deviations are that the in- 
volute strips form single curved surfaces and introduce a slight kink at the juncture of 
each surface-approximating element in the surface of revolution. These errors can be 
reduced to produce any degree of accuracy desired by increasing the number of surface- 
approximating elements and involute ribbon strips. 


DISCUSSION 

Applications to Heat -Exchanger Design 

There are many methods available for the approximate graphical development of 
double curved surfaces of revolution onto a plane (refs. 1 and 2). The method presented 
in this report differs from others in that a number of identical constant- width strips can 
be used to reproduce many surfaces of revolution. This feature should be of particular 
interest in the construction of heat exchangers where a imiform distribution of heat- 
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transfer tubes is required over a double curved surface of revolution. Two important 
structures of this type are tube -type rocket nozzles and the fuel- cooled nozzle plug for a 
jet engine. 

Regeneratively cooled tube-type rocket nozzles, where the involute ribbon strip 
method of development is used, can be constructed of constant- diameter tubing from the 
front end of the combustion chamber, through the throat, to the end of the expansion sec- 
tion. The number and size of tubes will be dictated by the requirements of critical heat- 
transfer regions in the nozzle assembly. The tube assembly, as in conventional designs, 
must be supported by a pressure-tight structure in order to contain the chamber pres- 
sure. 

Cooling a jet engine nozzle plug with liquid fuel can be accomplished by attaching 
cooling tubes to the shell which forms the plug. Since the plug is in a high-velocity, 
high-temperature gas stream, it is desirable to protect the coolant tubes, which gener- 
ally will be thin- walled and small in diameter. A direct approach is to install the tubes 
on the inside surface of the plug. One method for accomplishing this installation is to 
fasten the tubes to the developed involute ribbon strips for the plug shell before assembly 
into the shell form. In this way the fastening is simplified and inspection of the tube- 
strip joint is facilitated. If the tube- strip assembly is too difficult to shape over a man- 
drel, the tube and the strip can be preformed individually in a fixture to match a space 
curve identical to the composite involute on the surface of revolution. Assembly of the 
shell with the tube-strip elements would require a certain amount of springing over the 
mandrel for the plug. The final operation required to complete the plug would then be 
welding of the fa3dng surfaces of adjoining strips. 


Construction Problems 

Partially closed spherical shells and streamlined body shells or closed toroidal 
shells have a common assembly problem. In each case, some sort of internal removable 
support or mandrel is required for the retention of the involute strip dviring assembly 
and joining. The mandrel must be able to withstand welding or brazing operations. If a 
low-melting-point material is used for the mandrel structure, provisions must be made 
to insulate the material at those locations where welding or brazing is done. For a 
furnace-braze-type of joint, the mandrel could be made from materials such as plaster of 
paris which could be dissolved out of the assembly after joining. Since in all these struc- 
tures, with the exception of the annulus, two holes must be left in the shell (where the 
composite involutes start and end), it is conceivable, depending on hole size, that col- 
lapsible mandrels might be devised. In this event, the mandrels could be reusable. 

Where shell thicknesses are large, the involute strips would have to be preformed to 
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the final shape of the space curve, sprung into place over a temporary support, and tack 
welded to the adjacent strip. Before the shell is completely formed, the temporary sup- 
port may have to be removed, unless sufficient access was provided at the end holes. 


SUMMARY OF RESULTS 

The results of this investigation of the involute strip development method for fabrica- 
tion of rotationally S5unmetrical double cm-ved surfaces of revolution can be summarized 
as follows: 

1. Four types of rotationally S 3 nnmetrical surfaces of revolution have been developed 
by the involute ribbon method and presented in this report. These surfaces are (1) a par- 
tially closed streamlined body, (2) a rocket thrust chamber assembly, (3) a partially 
closed sphere, and (4) an annulus. The involute ribbons are the plane developed elements 
of the parent double curved surfaces of revolution. 

2. The generation of the specific composite involute ribbon for a double curved sur- 
face of revolution is dependent on the approximation of the surface shape by means of a 
series of contiguous cones, cylinders, and plane annulii, as appropriate. The greater 
the number of these approximating surface elements, the better the accuracy of the de- 
velopment. Since all these approximating elements are single curved surfaces, they can 
be developed exactly onto a plane surface. In a plane, one set of curves which are par- 
allel is a family of involutes. A developed helix is the limiting case for an involute. In 

a practical application, an involute ribbon would be used in place of the strip just defined. 
The composite involute curve obtained in going from one sirnface-approximating element 
to the next would then be the centerline of the ribbon. The development method can be 
accomplished graphically or by means of numerical computations. In the computational 
approach, the graphical method is reduced to the use of mathematical formulas where 
the composite involute curve is described in a suitable coordinate system. 

3. The computational approach produces better results than the graphical develop- 
ment method. In either method, accuracy can be improved by increasing the number of 
surface-approximating elements (cones and cylinders). 

4. The graphical involute development method can be very complicated because a 
multitude of required construction lines obfuscate the drawing. In addition, cones with 
small apex angles require the use of very long radii of curvature in the plane of develop- 
ment. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 16, 1970, 

720-03. 
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APPENDIX A 


DEVELOPMENT OF THEORY 
DERIVATION OF EQUATIONS 

The detailed derivation of the equations for developing a rotationally symmetrical 
surface of revolution is presented in the following text. The equations for defim'ng the 
centerline of a composite involute strip will be developed for a general surface of revo- 
lution which is composed of the basic surface elements of cone and cylinder. The plane 
annulus is a special case of the cone with a vertex angle of tt radians. In figure 8 are 
shown the basic elements for approximating any surface of revolution. By using many 
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Figure 8. - Basic elements for approximating any surface of revolution. 


small approximating elements, the actual surface of revolution can be approached as 
closely as desired. On the right side of figure 8, the cone (n + 1) is shown going either 
way so as to accommodate either an increasing or a decreasing slope. The left-hand 
cone, cone (n - 1), could have been shown this way too; however, the figure would then 
be too complicated. 


EQUATIONS FOR SINGLE CONICAL ELEMENT 

Before proceeding with the derivation of the equations for the general case, the 
method for a single conical element will be exposed. Figure 9 shows a developed ap- 
proximating surface in the form of a trimcated right- circular cone. This cone is labeled 
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Figure 9. - Developed conical surface and invoiute. 


n and is developed in a plane x-y. The x-y plane, representing rectangular Cartesian 
axes, is tangent to the cone surface, with the x-axis coincident with any arbitrary gen- 
erator of the cone surface. 

Coordinates of any point on the involute can be read directly from figure 9 as follows 


^ns ’?ns + ^ns 


(Al) 


^ns = ^ns ~ ^ns %sU 


or as 


^ns 

^ns ^ns^ 


(A2) 


For a particular width of involute strip, the surface area of the surface of revolu- 
tion can be obtained once the strip length is known. The length of the involute from the 
inner boundary of the developed cone and a point designated by S is obtained from fig- 

JZSo —2 

ure 9. A differential element of length is dS = t dx + dy and from equations (Al); 


= Pnd^ns ^^ns 
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whence 


S -- /”"" dins • ^ (lL - Ini) 

The area of the surface of revolution in figure 9 is very closely 

A = PioS 

Pno(’>ns-’>nl) 

A general surface of revolution, as approximated by the elements shown in figure 8, 
has been developed in figure 10, where for the sake of clarity the cone (n + 1) is assumed 
to decrease in diameter in going from left to right. In some cases, a plane annulus may 
be required as an approximating element for the surface, but this produces no complica- 
tions since this is a special case for a cone. The dimensions r and Au are, in gen- 
eral, readily available from a drawing or description of the surface of revolution. And 


rPn+1,0 



Figure 10. - Development of surface shown in figure 8. 
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the quantities p and cp are readily derived from r and Au as follows: 



(A5) 


and 


P 


ns 


ns 


sin (p 


n 


(A6) 


In this work the surfaces are rotationally symmetrical; therefore, the involutes will 
be generated from circles (the evolutes). The radii of the involute generating circles are 
determined from an arbitrary selection of strip width, the number of strips, and the re- 
quirement that a smooth transition occur in going from one surface-approximating ele- 
ment to its contiguous element. To preclude overlapping of the strip elements on the 
surface of revolution, it is necessary that the product v(ji not exceed the circumference 
of the smallest pole circle for the surface of revolution; that is, r^Q < r^^. 

The radius of the circle in the surface of revolution (see fig. 8) for determining the 
size of the involute generating circle for any cone n is 



(A7) 


Thus, r^^Q is a constant dependent only on the selection of v, oj and the requirement that 
adjacent involutes have a common tangent, which means that 


and 


^nl - '^n-1,2 
^n2 " ^n+1, 


(A8) 


where negative subscripts are excluded. The radius of the involute generating circle for 
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the n 


th 


surface element is 




n nO 

sin cp 


n 


^nO “ 


I/co 


2ir sin <p^ 


(A9) 


which is inversely proportional to sin <p^ for any particular selection of v and co. 

When is zero, p^q becomes infinite. For this case, = p^^ = , the surface 

element is a cylinder and the involute becomes a straight line in the plane of develop- 
ment. On the cylinder this line becomes a helix. For <p^ = n/2 radians, the surface 
element is a plane annulus. Therefore, in general, a surface of revolution can be de- 
veloped into a series of involute and helical strips. 

In figure 10 the composite involute curve is shown starting at the x-axis. The radius 
of the first involute generating circle is P^_^ q and is obtained from equation (A9), 
namely. 


p = I/Ct/ 

® 2)7 sin <p^ 1 

n-1 

The first segment of the composite involute, starting at the x-axis is extended through 
cone (n - 1). It is more convenient to use polar coordinates in defining the involute be- 
cause the radius vector to the involute also corresponds to a slant height along the cone. 
The definition of the angular argument y is readily obtained by reference to figure 9 and 
equation (A2). Equation (A8) implies that 

^nl = '^n-1,2 

~>'n2 " '^xv+1 , 1 


where negative subscripts are excluded. 

It will be shown that equation (A8), establishing the common tangency of contiguous 
elements of the composite involute curve, is a consequence of r^^Q being constant. By 
rearranging equation (A2) and exchanging j for s. 



(AlO) 
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and, therefore, 


n+(-l)j, 





n+ 


(- 1)^,0 


From figure 9, noting that j is exchanged for s. 


_ ^nj 


^nO ^nO 


and, similarly. 


^n+(-l)^, _ n+(-l)^, 

Pi r . 

n+(-l)%0 n+(- 1)^,0 

Also from figime 9 it can be seen that 

"’"'j ~'‘n+(-l)j,j-(-l)j 

and from equation (A7) 


Substituting equations (A14) and (A15) into equation (A13) yields 

^n+(-l)j, _ ^nj 

Equation (A16) substituted into equation (A12) shows that 


^n+(-l)j, _ ^nj 

^n+(-l)j,0 ^"0 


(All) 


(A12) 


(A13) 


(A14) 


(A15) 


(A16) 
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Thus, equations (AlO) and (All) are identical; and therefore, 

^nj =^n+(-l)j,j-(-l)0 


or 


^nl = Vl,2 

^n2 = '^n+1, 1 


(A8) 


EQUATIONS FOR SERIES OF CONICAL ELEMENTS OF 
MONOTONICALLY INCREASING SIZE 

The development, in the x-y plane, of a series of contiguous cones to the left of the 
cylinder in figure 8 is shown in figures 11(a) and (b). To be consistent with figure 8, the 
cones are labeled (n - 3), (n - 2), and (n - 1). Cone (n - 4) is the initial cone and is used 
here to establish a reference. Generally, the development can be started anywhere; 
however, it is convenient to start at an end cone (such as at the left-hand end). 


V 



(a) Cross section of part of surface of revolution for i itl2. 

Figure 11. - Development of series of contiguous cones to left of cylinder shown in figure 8. 


19 



■II 



(b) Composite involute generated on developed contiguous cones. 
Figure 11. - Concluded. 


From figure 11(b), it can be seen that successive values of the composite involute 
coordinates x and y can be written in terms of the immediately preceding values. 
Establishing the relations provides recursion formulas for determining points on the in- 
volute curves for each cone. From figure 11(b), it is apparent that 


^n-2,1 “Vs, 2’ V2,l “Vs, 2 


V2,2 “ Vl,l’ V2,2 “ Vl,l j 


(A17) 


Recalling that 


>^nl " Vl,2 
V " Vl, 1 
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and algebraically adding x and then y increments, the following equations are ob- 
tained: 


^n2 “ ^-1,2 " ^nl ^nl “ ^n2 ^n2^ 

^n2 = yn-1, 2 " <^nl ^nl " ^n2 ^n2^ 


(A18) 


EQUATIONS WHERE CONES ONLY ARE USED FOR APPROXIMATING 

A STREAMLINED SHAPE 

Before proceeding to the case where a cylinder is one of the surface -approximating 
elements, the case where the surface of revolution has a tear-drop shape and is approxi- 
mated by cones only will be considered. In the latter case, the composite involute curve 
will close (radii of curvature will decrease) when the approximating cones have a vertex 
half angle (p^ > n/2 because the radius vectors defining the involutes will cross over to 
the other side, as can be seen in figure 10 for cone (n + 1). 

In equation (A18), the algebraic sign of the second term is dependent on whether r^^g 
is greater or less than r^^j. Multiplying the second term of equations (A18) by 

^n l - ^n2 

I^nl-^n 2 l 

where the ratio 0/0 is defined to be zero, yields the equations for x^^g Yn2’ which 

for a tear-drop shape approximated by a series of colinear and contiguous cones become 

( ^nl " ^n2 \ 

T;: >^nl - Pn2 >^n2) 

l^nl " ^n2l' 

V (A19) 

^n2 = ynl + Vnl ^nl ' ^n2 >^n2) 

M^nl " ^n2l/ 
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EQUATIONS WHERE CONES AND CYLINDERS ARE USED 
FOR ANY SURFACE OF REVOLUTION 


The last step in the development of the equations for ^.nd the addition 

of a term to accoxint for a cylindrical element in the surface of revolution. The incre- 
ments in X and y can be read directly from the diagram of figure 12. The developed 


y 



Figure 12. - Development of cone-cylinder-cone surface of revolution showing composite involute curve. 


heUx is tangent to the connecting involutes at each end, establishing its slope in the x-y 
planes as the tangent of the angle From the figure then 

AXjj = 5 cos 

sin^<ic^ 


Ay = ^ sin (aci 
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The magnitudes of angle and angle ao<f can be read directly from figure 12 as 


= JT - 

(^ac<f -l^ck - {aclt’ 


^nl + ^nl 


In the determination of the x and y increments for the cylindrical element, 


^nl = ^n2 


and 


>^nl = ^n2 


whence, 




-h„ cos 
cos(T,ni - rni) 


^n ^nl 

cos(T,ni - y^i) 


The denominators of these equations can be modified by the use of equation (A2); then 


^^n = 


Ay = 

•^n 


-hn cos ^nl 

cos arctan 

^n ^nl 
cos arctan 


(A20) 


The addition of the x and y increments from equation (A20) to eq\iation (A19) will give 
the values of x and y when a cylinder also is used in the approximation of a surface of 
revolution. When this addition is made, the generality of the equation can be preserved 
by multiplying the term contributed by the cylinder by an operator f^ which has the 
following values: 
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0 

II 

when 

■"ill * ■•n2 

II 

when 



Making these substitutions results in a pair of general recursive equations for x and y 
for developing double curved surfaces of revolution, namely, 


^n2 




/ ^nl - ^n2 \ 

\l^nl ■ ^n2l>' 


^nl ’^nl - Pn2 ^n2^ + ^n 


\ ^nl 




cos arctan 


^nl 


> (A21) 


^n2 = ynl + ( -T )^nl y„i - P„2 sin 7 ^ 3 ) + 

\ ^nl “ ^n2 ' 


hn sin 


cos arctan 


^nl 


For the determination of points on the composite involute curve at locations, on the cylin- 
der, designated by the subscript ns, the equations are modified to read as follows: 


X, 


ns 


'ns 


^1 


^nl 


^nl - ^n2 
^nl - ^n2 


Pnl ^nl - ^n2 7^2^ + ^n 


/ ^nl - ^n2 \ 

Vl^nl-^n2l/ 


Kl ^nl - ^n 2 ^n 2 ^ + ^ 


^%s ^nl 
cos arctan 

^nl 


> (A22) 


Au^s sin 


n 


cos arctan 


^nl 


LENGTH OF COMPOSITE INVOLUTE CURVE 


A formula for computing the length of the composite involute curve to any location 
denoted by the subscript ns is developed here. Equation (A3) can be modified by the 
substitution suggested by equation (A2), where 



then 




(A23) 
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For a set of contiguous involutes plus a segment of the n 


th 


involute 




1 (pIs - pIi \ 
A ^nO / 


i=l 


(A24) 


Complications arise where a helix, resulting from the use of a cylinder as an approxi- 
mating element, is interspersed in the composite curve. Referring to figure 12 it is 
seen that the length of the helical element is 


cos arctan 

Since it may be desired to obtain the length of the helix to any point designated by the 
subscript s on the helix, the segment of the cylinder is more conveniently defined by 
substituting for h. Thus the length of helix for the n^" cylinder is 

u - u 1 
S - — nl 

^ cos arctan ri^-^ 


In order to combine the lengths of the helical elements with the lengths of the involute 
elements, use of the operators applied in equations (A22) can be made here, namely. 


^nl - ^n2 

l^nl - ^n2l 


and f 

n 


0 

1 


when r^i ^ r^2 
when r^i = r^2 


For application to equation (A24), 


the first operator can be written as 


^il ~ ^i2 

l^il - ^i2 


and 


^ns - ^n2 
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where 0/0 is defined as equal to zero. Tlie second operator can be rewritten as 
where 



when r .2 

when r. ^ = r -2 


These operators have the purpose, as in equations (A21) and (A22) of admitting only those 
terms which represent actual length increments. 

Combining the lengths of all the involute elements with the lengths of all the helical 
elements yields the following formula for the length of the composite curve to any de- 
sired location: 



(A25) 


noting that 77n+i, 1 = % 2 ‘ 

TRACE OF COMPOSITE INVOLUTE CURVE ON SURFACE OF REVOLUTION 

Formula.s defining the trace of the composite involute curve on the surface of revo- 
lution are developed here. The axes are labeled u, v, and w in cyclic order for a 
right-hand orthogonal system. The u-axis is made coincident with the axis of S}rmmetry 
of the surface of revolution. These axes and the other pertinent dimensions are shown 
in figure 13. Figure 13(b) shows the first two approximating cone elements developed in 
the x-y plane, as defined previously, where the x-axis is located in the u-w plane and the 
x-y origin is located at the beginning of the first conical element. For the sake of clar- 
ity, the addition of a cylindrical element will be made later. 

The u, V, w coordinates are obtained as functions of p, <p , and 9. In order to 
simplify the derivation, the value of is assumed to be equal to zero. This gives 
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Cone 2 


(a) Two contiguous cones approximating a surface of revolution. Curve kac is involute 
for cone 1; curve nfcd is involute for cone 2. 



(b) Planar development of surface shown in figure 13(a). Corresponding points have the same notation. 
Figure 13. - Development of composite involute for general surface of revolution approximated by cones. 


the trace for one composite curve which starts with the coordinates 


^11 "^11 


Wii sin cp^ 

The traces of other composite involutes which would be the centerlines of other strips 
can be obtained by setting = 27rn / v radians, the initial angle about the u-axis for the 
n^^ strip. 

In figure 13, corresponding points on the surface of revolution and on the plane of 
development are labeled identically. Equivalence of arc lengths in the u, v, w coordi- 
nate system and the x-y plane will be used to calculate the values of 0 and its incre- 
ments ij/. From figure 13(a) 

^ C - Pi2^12 ^1 

and from figure 13(b) 

io =Pi2(ri2 - rji) 


whence 


^12 - ni 

sin (pj 


For continuity of values of 0 


®n+l, 1 - ^n2 


Similarly, for the arc di’ from figure 13(a) 

^ ^ 22^22 ^2 


and from figure 13 (b) 

= ^ 22^^22 “ '^ 21 ^ 
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•whence 


Note that since = 0, 
that 


and further that 


and that, generally, 


or written in terms of 


and 


for intermediate points. 


'>'22 ~ "^21 
sin <p 2 


then = ^^21 ^ ®12- 


Again from figure 13(a), it can be seen 


^22 


'>'22 ~ ’>'21 
sin ^2 


>'12 - >'11 >'22 - >'21 

— — — — + 

sin cp^ sin cp 2 


^32 “ ^22 *^32 

>'12 ~ >'11 ^ >'22 ~ >'21 ^ >'32 ~ >'31 
sin cp ^ sin cp 2 sin q) ^ 


"n2 


= Vl,2 + 


n2 


and (p 



n-2+j 

E >'i2 - >'il 

sin 


Z >'i2 ~ >'il ^ >'ns ~ >'nl 
sin cp^ sin cp^ 
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Figure 14. - Surface of revolution in form of right -circular cylinder. 


A formula for 9^2 written in more general form where, in addition to 

the cones used in the approximation, cylinders can also be used. In figure 14 a simface- 
approximating element in the form of a right- circular cylinder is shown. The same 
cylinder in the plane of development is shown in figure 12 as element n. From figure 14 



’^n2 


and from figure 12 



h 


n 


cos(77ni - r„i) 



whence 

fd = h^ tan(r?^l - 


From equation (A2) it is seen that 


%1 - ’'nl = ’'nl 


therefore 


fd = h^ tan arctan 
= \ ^nl 
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The arc (frd is equal to the straight-line segment then 

’^n2^n2 = ^n^nl 


whence 


h 




n 


n2 


n2 


^nl 


This is the increment in the argument 6 in going from one end of the cylinder, along a 
helix, to the other. 

An equation for 6^2 which includes both conical and cylindrical elements can now 
be written. The operators used in equations (A21) can be used to advantage here, re- 
sulting in 


'n2 


i=l 


/ -'ll - >-12 \ 

M 2 - -yn) 

bf ^ 

Vl 

^il-^i2l^ 

\ sin ) 



(A26) 


Referring to figure 13(a), equations for the u, v, w coordinates of the composite 
involute can now be written as the following set of recursive relations. 


“n2 * “nl + ( ‘ «’n " 

^ ^nl " ^n2l' 


^n2 = ^n2 ^n2 


^^n2 = ^n2 ^n2 


(A27) 


J 
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APPENDIX B 


COMPUTER PROGRAM 

A geometry calculation code, STRIP, was programed to compute the planar coordi- 
nates needed to define the basic involute ribbon of width w. N number of these ribbons 
are used to fabricate a particular rotationally symmetric surface by a method described 
in this report. Specific equations used by the program are shown in the order of their 
use in the program. Specifications of input data are given. Outputs from a sample prob- 
lem are included to demonstrate the program usage . The problem was programed en- 
tirely in FORTRAN IV language for a 7044/7094 DCS computer. 


PROBLEM DESCRIPTION 

To solve for the plane coordinates that define a ribbon, the body must be divided into 
a number of cones along the axis of symmetry. Each cone is numbered consecutively 
from left to right and is associated with a left cone radius and the cone height (see 
fig. 15). The cylinder is considered to be a special case of the cone. When a large cone 




height is specified, up to eight intermediate cones can be specified to help smooth the 
ribbon profile. 


Equations 

The calculations in this program concern themselves with defining the ribbon in 
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terms of p, y, x, and y. The equations used are developed in appendix A and rear- 
ranged in this appendix to facilitate programing. Separate equations for the cylinder 
case were developed. 

Cone equations . - The cone equations are as follows: 

(1) Cone constants: 


(2) Polar data: 


(3) Cartesian data: 


sin cp„ = 


lyi-^n+l,ll 




nO 


U(jO 


2n sin co 
^ n 


n, m 


"n,l 

sin tp 


n 




y = T) - arctan rt 
^n, m 'n,m 


n, m 


^n, m ^n-1, m ^ ^n, 1 ^n, 1 ^n, m ^n, m 

^n, m “ ^n-1, m ^ ^n, 1 ^n, 1 ^n, m ^n, m 


(Bl) 


(B2) 


(B3) 


(B4) 

(B5) 

(B6) 

(B7) 


where m designates the intermediate radii for cone n. No intermediate radii are 
allowed for the cylinder. 

C ylinder equations . - The cylinder equations are as follows: 

(1) Polar data, for the cylinder in the initial position (n = 1): 


77 1 1 = arcsin/ — — — ^ 

V2.r,i/ 

n,l =r,i i - arctan rji^. 


(B8) 

(B9) 
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(2) Polar data, for the cylinder in any position except initial position: 


Vl - V2 


^n, 1 ^n-l,m ^n, 2 


(3) Cartesian data, for the cylinder in the initial position: 

xi,i=0 


h ,2 -^ 1,2 ^ 1 , 2 , 


yi,i 


= 0 


^ 1,2 -^ 1,2 ^ 1 , 2 ^ 

(4) Cartesian data, for the cylinder in any position: 

h_ 


V2 


n 


®‘"(f + ^'n, 1 - ”n, 1 
>'n,2 


X. 




(BIO) 

(BID 


(B12) 


(B13) 


(B14) 

(B15) 

(B16) 


Input Data 

The input into STRIP consists of three card formats. The information required to 
enter into these cards is described here. Included in the description is the number of 
input cards required. The actual formats of the cards follow the description. 

Description. - The following information is required to enter into the STRIP card 
formats: 

(1) Identification: The first card of each case contains a descriptive heading which 
will appear in the output as a heading in various places. 

(2) Ribbon data: The second card specifies 

(a) The total number of radii starting from left to right associated with each 

cone or cylinder, including the final closing radius 

(b) The total number of ribbons 
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(c) The width of ribbon (If the given width is too large to be consistent with the 
smallest radius, the program will compute and use a corrected ribbon 
width. ) 

(3) Cone data: The third and succeeding cards describe the model as it is simulated 
by contiguous cones. For each cone or cylinder, the following data are specified. 

(a) The cone index number 

(b) The number of intermediate equal-height cones to be added (A maximum 

number of eight cones may be added. For the cylinder (a special case of 
the cone) case, a 1 is entered. For the final most right-hand radius, a 0 
(zero) is entered. ) 

(c) The actual value of the radius 

(d) The cone height 

Data for two cones may be specified on each card. A card need not be full but the 
first set of columns must be filled. The read operation ends with a blank card or a 
card with the first data field blank. Cases may be stacked. Each subsequent case 
starts with an identification card. The input data are independent of the system of length 
units used. 

Card formats . - Examples of card formats are 

(1) Card 1: Identification (Format 9A6) 

1 

54 columns may be used for a descriptive heading 

(2) Card 2: Ribbon data (Format 215, FIO.O) 


1 5 

6 10 

11 20 


Number of 

Number of 

Width of ribbon, w 


radii, NR 

ribbons, v 



X 

X 

x.xxxx 



(3) Card 3: Cone radius and heat data (Format 2 (215, 2F10. 0)) 


5 

10 

20 

30 

35 

40 

50 

60 

Cone index, 
n 

Number of 
intermediate 
radii, m 

Cone radius, 
^n,l 

Cone height, 

hn 

n 

m 

Vi 


X 

X 

X.XXXX 

X.XXXX 

X 

X 

X.XXXX 

X.XXXX 
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(4) Blank card (last card) 


In these formats X represents a numerical digit: 

(1) X indicates that the field contains an integer and is right justified in the card 
field. 

(2) X. XXXX indicates a fixed decimal point number with the decimal point occur- 
ring anywhere in the card field except the first position which is usually reserved for 
the sign of the number. 

Problem lim its. - There are size restrictions which must be considered due to 
core storage size. The maximum number of cones is set at 100. Due to printing limi- 
tations, eight intermediate cone radii may be specified to smooth ribbon data for cones 
with large cone heights. 

Cone number ing. - Numbering starts with index number 1 and must be sequential 
to the maximum number 100. To facilitate changing the model, the cone data cards can 
be in any order as long as each card has an index number in the first card field. H two 
of more cone data cards have the same index, the last cone data card read will be used. 


Output Data 

The printed output consists of the following: 

(1) The input 

(2) Complete cone data including intermediate radii 

(3) Computed polar data, polar radius and angle 

(4) Computed x-y coordinate data of the ribbon centerline 

(5) Strip length against model length 


EXAMPLE 

A listing of the program with one complete example of a jet engine nozzle plug 
(figs. 3 and 16) is presented which includes a listing of input cards and a computer listing 
of the output data. 
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Figure 16. - Input data for sample problem (jet engine nozzle plug). 


Computer Program STRIP: Program Listing 


$IBFTC STRIP 


C 

C 

C 

C 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


«**INPLT CAROS INSTRUCTIHNS+'fr* 

CARO ). - TITLE CARO 

FORMAT 9A6 (54 ACTIVE COLUMNS) 

CARD 2 - RIRfiON DATA 

FORMAT 21 5,F10.0 

NR = NO, OF RAOI I 

NL = NO. OF RIBRONS(STRIPS) 

W = WIDTH OF RIBBON (CORRECTED BY PROGRAM IE TOO LARGE) 
CARO 3 TO NR + 3 - CONf RADIUS AND HEIP.HF DATA 

FORMAT 2 ( 21 5, 2F10.0) (FIRgT F I El D BLANK OR BLANK CARD ENDS 

READ ). OR 2 CONE DATA OLR CARD) 

I = COVE INDEX NO. FROM L TO R 

J = NO. OF INTERMEDIATE RADII TO BE ADDED, MAX OF 8 
(J=l FOR CYLINDER, J=0 FOR LAST RADIUS) 

R = LEFT RADIUS OF CONE 
H = CONE HEIGHT 

-REPEAT ABOVE FOR SECOND SET ON CARD- 
CASES MAY BE STACKED. NEXT CASE STARTS WITH TITLE CARD. 
COMMON/IDENT /I0(9) 

COMMON/ INPUT /R(lOi,9),H(100),M( 100 ) ,NU, W, NC 
COMMCN/MIOATA/P ( 1 00 ) , RO ( 1 00 , 10 ) , E T A ( 1 00 , 10 ) , G AM ( 100 , 10 ) 

COMMON/ OUT PUT /X( 100, 10), Y( 100, 10) 

REAL MINR 


37 



C INITIALIZATIONS 
C 

300 NMAX = 0 

DO 1 1=1,100 
DO 2 J=l,10 
R( I , J ) = 0.0 
R0( I , J) = 0.0 
ETA ( I ,J ) = 0.0 
GAM( I , J ) = 0.0 
X( I , J) = 0.0 
Y( I,J) = 0.0 
2 CONTINUE 
H( I ) = 0.0 
M( I ) = 0.0 
P( I ) = 0.0 
1 CONTINUE 

TPI = 2.0*3.1415926 
C 

C INPUT/OUTPUT ID AND PROBLEM DATA 
C 

READ(5,901) ID 
WRITE (6 ,900) 

WRITE{6,902) ID 
READ(5,903) NC,NU,W 
WRITE (6, 905) NC , NU , W 

10 READ ( 5, 904 ) I , J , A , B , K , L , C , D 
IF( l.EQ.O) GO TO 11 

NMAX = MAXOINMAX , I ,K) 

M( I ) = J 
R( I , 1 ) =A 
H( I ) = B 

IF(K.EQ.O) GO ID 10 

M(K) = L 

R(K, 1 ) = C 

H(K) = 0 

GO TO 10 

11 CONTINUE 
C 

C CHECK FOR CONSISTENCY 
C 

RU = NU 

MINR = W*RU/TPI 
CHECKR = R(l,l) 

DO 5 1=1, NMAX 

CHECKR = AMIN1(CHECKR,R( 1,1) ) 

5 CONTINUE 

IFIMINR.LT. CHECKR) GO TO 6 
W = CHECKR*TP I/RU-.OOl 
WRITE (6,915) W 

6 NMAX = NMAX-1 
C 

C COMPUTE COMPLETE R MATRIX 
C 

DO 12 1=1, NMAX 
L = M( I ) 

R I = R ( I , 1 ) 

RJ = R( 1+1,1) 

RJMI = RJ-RI 
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DEBUG KI ,RJ,RJMI ,L 
00 13 J=1,L 
A = J 
B = L 

SLOPE = A/B 

R(I,J+1) = RI+SL0PE*RJMI 
13 CONTINUE 
12 CONTINUE 

OUTPUT COMPUTED PLUS INPUT DATA 

WRITE(f>,906) 

DO lA I=1,NMAX 
L = M( I ) + l 

WRITC(5,913) I,R(I,l),R(I + l,l),Hm,(R(I,J),J = 2,L) 

lA CONTINUE 

COMPUTE POLAR DATA 

DO IS I=1,NMAX 
L = M( I )+l 

RJMI = ABS(R( H-1 , 1)-R( I , 1 ) ) 

IFIRJMI .EQ.0.0) GO TO 17 

SPHI = RJMI/SCJRT (H( I )<'’f‘2 + RJMH<«2) 

RNO = RU^i-W/ ( TPI «SPOI ) 

DO 16 J=1,L 

R0{ I , J) = R ( I , J ) /SPHI 

ETAII.J) = S(yRT( (RU( I,J)/RN0)f'«2-1.0) 

GAM(I,J) = ETAd.JI- ATAN (ETA(I,J)» 

16 CONTINUE 
GO TO 15 

CYLINDER CASE 

17 L = M(I-1)+1 

R0( 1 , 1) = 0.0 

IF(I.EQ.l) GO TO 18 
ETA ( 1,1) = ETA( I -1 ,L ) 

GAM( 1,1) = GAM( I -1 , L ) 

GO TO 19 

18 £TA(I,1) = ARSIN (RO*W/( TPI*R ( I , 1 ) ) ) 

GAM(1,1) = ETA( I , 1 I-ATAN ( ETA ( I , 1 ) ) 

19 RO( I ,2) = R0( 1 ,1 ) 

ETA{ 1,2) = ETA( 1,1) 

GAM( 1,2) = GAM( 1,1) 

15 CONTINUE 

COMPUTE X-Y COORDINATE DATA 

DO 20 I=1,NMAX 
L = M ( I ) 

TEST = R(I+1,1)-R(!,1) 

IF(TEST) 21,22,23 

21 Q = 1.0 
F = 0.0 
GO Tc: 2 A 

22 Q - 0.0 
F = 1.0 
GO TO 2A 
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23 Q = -1.0 
F = C.O 

24 RONl = R0( If 1) 

RCG = RGN1*C0S( GAM(I , 1) ) 

RSG = RQN1*SIN(3AM( I , 1) ) 

IF(I.NE.l) GO TO 25 
XCIfl) = RCG 
Y( 1 , 1) = RSG 
GO rC 26 

25 K = M(I-1)+1 

X( I , 1 ) = X( l-l,K ) 

Y<I,1) = Y(I-1,K) 

26 DO 30 J=1,L 
IFIF.EQ. 1.0) GO TO 27 
0 = R0( 1 , J + 1) 

A = GAM( I , J+1) 

DX = ( RCG-0*C0S ( A ) ) 

OY = Q<'(RSG-D#SIN(A) ) 

GO TO 28 

27 D = ABS (H( I ) /SI N ( TPl /4. 0+GAM ( I , J )-ETA( I , J ) ) ) 

A = TPI/2.0-ETA( IfJ) 

OX = -D*COS(A) 

IF( I .EO . 1 ) DX = -DX 
DY = D*S1N(A) 

28 X ( I , J + 1 ) = X ( I, 1 )+DX 
Y( I , J+1 ) = Y( I, 1 )+DY 

30 CONTINUE 
20 CONTINUE 

OUTPUT COMPUTED POLAR DATA AND X-Y COORDINATE DATA 

WRITE (6, 907) 

WRITE(6,908) 

DO ICO I=1,NMAX 
L = Mil) 

WR ITE (6 , 914) I ,R0( I , 1 ) ,R0( I , L + 1 ) , ( R0( I , J ) , J = 2,L ) 

100 CUNTINUE 
WRITE (6 ,909) 

WRITE(6,908) 

DO 101 I=1,NMAX 
L = Mil) 

WRITC(6,914) IfGAMIIfl) ,GAM( I ,L + 1 ) , (GAM( I , J ) , J = 2,L ) 

101 CONTINUE 
WRITE (6,910) 

WRITEI6 ,908) 

DO 102 I=1,NMAX 
L = M( I ) 

WRITE(6,914) I,X(I,1),X(I,L+1),{X(I,J) , J = 2 ,L) 

102 CUNTINUE 
WRITF(6,912) 

WRITE(6,908) 

DO 103 I=1,NMAX 
L = M(I ) 

WRITE(6,914) I,Y(I,1),Y(I,L+1),(Y(I,J),J=2,L) 

103 CONTINUE 
GO TO 300 

900 FORMAT! IHl, 5X,85HSTRIP - A PROGRAM TO APPROXIMATE A BODY OF REVO 
ILUTION BY A SERIES OF CONTIGUOUS CONES/ 

2 5X, 36H0GVELGPE0 FROM INVOLUTE CURVE STRIPS) 
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901 F0RMAT(9A6) 

902 FORHATdHO, 5X,1AA6) 

903 F0RMAT(2I5 ,F10. 0,215) 

904 F0RMAT(2(2I5,2F10.0) ) 

905 FORMAT! IHO, 5X,16HFIXED INPUT DATA/ 

1 5X,18HNUMB£R 3F CONES = ,15/ 

2 5X,39HNUMBER 3F STRIPS ON THE CIRCUMFERENCE = ,15/ 

3 5X,29HRIBB0N WIDTH OF EACH STRIP = ,F10.5) 

906 FORMATdHO, 5X , 1 9HVAR I A B LE INPUT DATA/ 

11X,4HC'JNE,6X, lOHCONE RADII, 8X, 4HC0NE , 30X , 18H I NTERMEO I ATE RADII/ 
22X,2HN3, 6X, H 1 ,9X ,192 , 9X,2HHT, 8X , IH 1 , 9X , 1H2 , 9X , 1H3 , 9X , 1 H4 , 9X , 

31H5,9X,iH6,9X,lH7,9X,lH8//) 

907 FORMAT! IHO, 5X, 1 9HCUMPUTED POLAR DATA/ 

11X,4HC3NE,6X, IIHPOLAR RADII, 7X, 34X , 16H I N T ERMED I ATE RHO) 

90 8 FORMAT! 2 X , 2HN0, 6 X , 1 HI , 9 X , 1H2 , 19X , IHl , 9X , 1H2 , 9X , 1H3,9X, 1H4,9X, 1H5, 
19X,1H6,9X, lH7,9X,lfi8//) 

909 FORMAT! IHO, 5 X , 1 9HC 0 MPU T ED POLAR DATA/ 

11X,4HC3NE,6X,11HP0LAR ANGLE, 7X, 34X , 18H I NT ERMED I ATE GAMMA) 

910 FORMAT! IHO, 5X , 2 3HC OMPU T ED X-Y COORDINATE DATA/ 

11X,4HC3NE , 6X, 12 9 /-COORDINATE , 6X , 34X , 19H I NT ERMED I ATE X-OATA) 

912 FORMAT! IHO, 5 X , 2 RHC OMPU T ED X-Y COORDINATE DATA/ 

11X,4HC JNE ,6X , 129 Y-CUORDI NATE ,6X , 34X , 19H I NT ERMED I ATE Y-DATA ) 

913 FORMAT! 1H0,I4,1X,11F10.5) 

914 FORMAT 1 IHO , I 4 , IX , 2F 10.5 , lOH *4<<=****>!= , 6 F 1 0 . 5 ) 

915 FORMAT! 5X,39HC3RRECTED RIBBON WIDTH UF EACH STRIP = ,F10.5) 

END 


FORTRAN 

symbol 

ETA 

GAM 

H 

ID 

NR 

NU 

R 

RO 

W 

X 

Y 


STRIP Dictionary 

Definition or equivalent 
mathematical symbol 

V 

ribbon polar coordinate, y 

cone cylinder height, h 

title and problem identification data 

total number of radii from model data 

u 

model radii, r 

ribbon polar coordinate, p 

cu 

ribbon Cartesian coordinate, x 
ribbon Cartesian coordinate, y 



Listing of Input Cards 


SAMPLE PROBLEM - 
13 50 .25 

I b~2 .V9 

3 8 A. 17620 


7 8 6.5A861 

“5 TTVrSVSI 

11 8 7.73482 

T3 8"" 5 '.'77 


PLUG 

3.142T' 

3.1425 

“3il425 

3.1425 

1.23 

2.0 

TT49 


2 ■ 8 3.58310 

4 8 4.76931 

6 6 5.95551 

8 8 7.14172 

'“10 1 7.73482 

12 8 7.0 

14 2.14 


'3.1425 

3.1425 

3T14Z5 

3.1425 

1.23 

1.59 


Listing of Output Data 


STRIP - A PROGRAM m APPROXIMATE A BODY OF REVOLUTION BY A SERIFS UF CONTIGUOUS CONES 
DEVELOPEJ FROM INVOLUTE CURVE STRIPS 

SAMPLE PROBLEM - PLUG 

FIXED INPUT DATA 
NUMBER OF CONES = 13 

NUMBER OF STRIPS ON THE CIRCUMFERENCE = 50 

RIBBON WIDTH OF EACH STRIP = J. 25000 

VARIABLE INPUT DATA 


CONE 

CONE 

RADI I 

CONE 



INTERMEDIATE 

RADI I 




NO 

1 

2 

HT 

1 ' 

2 

3 

A 

5 

6 

7 

8 

L 

2.99U00 

3.58310 

3. 1A250 

3-06A1A 

3.13827 

3.212AI 

3.28655 

3.36069 

3.A3A82 

3.50895 

3.58310 

2 

3.58310 

A. 17620 

3. 1A250 

3.6572A 

3-73137 

3.80551 

3.87965 

3-95379 

A. 02792 

A. 10205 

A- 17620 

3 

A. 17620 

A .76931 

3. 1A250 

A.2503A 

A-32AA8 

A. 39862 

A.A7275 

A-5A689 

A. 62103 

A. 69517 

^1, 76931 

A 

A. 76931 

5.362A1 

3.1A250 

A.0A3A5 

A. 91756 

A-99172 

5.06586 

5. lAOOO 

5.21A13 

5.28827 

5.362A1 

5 

5.362AI 

5.95551 

3. 1A250 

5.A3655 

5.51068 

5-58A82 

5.65896 

5. 73310 

5.8072A 

5.88137 

5.95551 

6 

5.95551 

6.5A861 

3 . 1A25Q 

6.02965 

6-10378 

6.17792 

6.25206 

6.32620 

6.A0O33 

6. A7AA7 

6.5AB61 

7 

6 • 5 A8 6 1 

7.1A172 

3 . IA250 

6.62275 

6.69689 

6.77103 

6.8A516 

6. U930 

6.993AA 

7.06753 

7.1A172 

8 

7. lAl 72 

7.73A82 

3 . 1A250 

7.21586 

7-28999 

7.36A13 

7. A3827 

7.N12A1 

7.5865A 

7.66068 

7.73A82 

9 

7. 73A82 

7.73A82 

1 .23000 

7.73A82 








10 

7. 73A3? 

7.73A32 

1.23000 

7.73A82 








11 

7.73A82 

7.0000U 

2.00000 

7.6A297 

7.55111 

7.A5926 

7.367A1 

7.27556 

7. 18370 

7.09185 

7.00000 

12 

7. CQOOO 

5.77000 

1.59000 

6.8A625 

6.69250 

6-53875 

6-38500 

6.23125 

6.07750 

5.92375 

5.77000 

13 

5. 7 7000 

2.1A000 

1. A9000 

5.31625 

A-86250 

A.A0875 

3.95500 

3. ‘>0125 

3.0A750 

2.59375 

2 . lAOOO 


COMPUTED POLAR DATA 










CONE 

POLAR 

RADII 




INTERMEDIATE RHU 




NO 

1 

2 


1 

2 

3 

A 

5 

6 

7 

8 

1 

16. 12200 

19.3199b 

# ^ « « 4 = ^ 

16.52175 

16.92150 

17-3212A 

17.72099 

18. 1207A 

18.520A9 

18.92023 


2 

19. 31998 

22.51 796 


19.71973 

20. 119A8 

20.51922 

20.91897 

21. J1872 

21. 71HA7 

22.11821 


3 

22.51759 

25.71557 


22.9173A 

23-31709 

23.7168A 

2A, 11658 

2A. 51633 

2A. 91608 

25.31583 


A 

25. 71590 

28.91397 


26. 1 157A 

26.515A9 

26.9152A 

27.31A98 

27. 71A73 

28.ll'iA8 

28.51A23 


5 

28.91397 

32.11195 


29.31372 

29.713A7 

30.11322 

30.51296 

30.91271 

31.312A6 

31.71221 


6 

32.11196 

35.3099A 


32-51170 

32.911A5 

33.31120 

33. 71U95 

3A. 11069 

3A.510AA 

3A. 91019 


7 

35. 30936 

38.5073A 

4 « * 4c 4 4 4 * 

35.70910 

36. 10885 

36.50860 

36.9D835 

37.30810 

37.7078A 

38.10759 
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8 

38.50797 

A1 .70595 


38.90771 

39.307A6 

39.70721 

AO. 10696 

AO- 50670 

A0.906A5 

Al. 30623 

9 

0. 

0. 


0. 







10 

0. 

0. 


0. 







11 

22.A282A 

20-29753 


22. 16190 

21.89556 

21.62922 

21.36288 

21.0965A 

20-83021 

20.56387 

12 

11.AA030 

9.A3007 


11.18902 

10.9377A 

10.6B6A6 

10.A3519 

10. 18391 

9.93263 

9.68135 

13 

6.23717 

2.31 326 


5.7A668 

5.25619 

A. 76570 

A. 27521 

3. 78A73 

3.29A2A 

2.00375 

CONE 

NO 

COMPUTED POLAR DATA 
POLAR ANGLE 
1 2 


1 

2 

INTERMEDIATE GAMMA 
3 A 5 

6 

7 

1 

0.27916 

0.51578 


0.3072A 

0.33583 

0.36A87 

0.39A33 

0.A2A19 

0-A5A39 

0.A8A93 

2 

0.51578 

0.771A2 


0.5A691 

0.57831 

0.60996 

0.6A18A 

0.6739A 

0.70625 

0.73875 

3 

0.771A2 

1.03828 


0.80A27 

0.83729 

0.870A5 

0.90376 

0.93720 

0.97077 

l.OOAAT 

A 

1.C3828 

1-31236 


1.07221 

1.10623 

1.1A036 

1. 17A59 

1.20891 

1.2A331 

1.27783 

5 

1. 31236 

1.591A0 


1.3A701 

1.38173 

1.A1651 

1.A5137 

1.A8628 

1.52127 

1.55631 

6 

1.591A0 

1.87399 


1.62655 

1.66176 

1.69701 

1.73232 

1.76767 

1.80307 

1.83851 

7 

1.87399 

2-15923 


1.90952 

1.9A508 

1.98069 

2.01633 

2-05200 

2.08771 

2- 123A5 

8 

2.15923 

2,<»A6A7 


2. 19503 

2.23087 

2.2667A 

2. 30263 

2. 33855 

2.37A50 

2.A10A7 

9 

2. AA6A7 

2,«»A6A7 


2.AA6A 7 







10 

2.AA6A7 

2.AA6A7 

*«r*4>**«* 

2. AA6A7 







11 

2.AA6A7 

2.09087 


2.A0188 

2.35732 

2.31280 

2.26832 

2.22389 

2.17950 

2.13516 

12 

2.090R7 

1.50369 

«««««««« 

2.01685 

1.9A298 

1.86927 

1.7957A 

1. /22AD 

U6A927 

1.57636 

13 

1.50369 

O.Ol 899 

««««*««« 

1.29083 

1.0809A 

0.87A99 

0.67AA7 

0.A817A 

0.30089 

0.13991 

CONE 

NO 

COMPUTED X-Y COORDINATE DATA 
X-COOROINATE 
1 2 

1 

2 

INTERMEDIATE : 
3 A 

X-DATA 

5 

6 

7 

L 

15.A9789 

16.80660 


15.7A806 

15.97623 

16. 18099 

16.36095 

16.51A7B 

16.6A116 

16.738BA 

2 

16.80660 

16.1A35A 


16-8A32fl 

16.8A775 

16.81898 

16.7559A 

16.65772 

16.523AA 

16.35229 

3 

16. 1A35A 

13.05608 


15.89651 

15.61063 

15.28537 

lA. 92031 

lA. 51508 

1A.069A1 

13.58311 

A 

13.05608 

7.389AA 


12-A8829 

11.87982 

11,23081 

10.5A151 

9.81226 

9. 0A3A7 

8-23567 

5 

7. 389AA 

-0.66160 


6.505A9 

5.58A60 

A.6276A 

3.63556 

2 .609A2 

1.55036 

0.A5953 

6 

-0,66160 

-10.5A255 


-1.81179 

-2.98953 

-A, 19328 

-5.A21A0 

-6. -7219 

-7.9A389 

-9.23A6A 

7 

-10.5A255 

-21.37392 


-11.86562 

-13.2018A 

-IA.5A909 

-15.90523 

-17.26806 

-18.63532 

-20.00A72 

8 

-21.37392 

-32.02901 


-22.7A053 

-2A. 10215 

-25.A5635 

-26.80065 

-28.13253 

-29.AA963 

-30.7A92B 

9 

-32.02901 

-35.933A5 


-35.933A5 







10 

-35.933A5 

-39.83789 


-39.83789 







11 

-39.83789 

-A6.9755A 


-AO. 69217 

-A1.562A0 

-A2.AA600 

-A3.3A039 

-AA.2A305 

-A5.151AA 

-A6. 06308 

12 

-A6.9755A 

-53.29308 


-A7. 83372 

-A8.6832A 

-A9. 51823 

-50.33311 

-51.12267 

-51.88207 

-52.60687 

13 

-53.29308 

-55.18768 

*#««**** 

-5A.A6276 

-55.3A787 

-55.92966 

-56.2139A 

-56.22882 

-56.02107 

-55.65119 


COMPUTED X-Y COORDINATE DATA 


CONE 

Y-CQOROINATE 




INTERMEDIATE Y 

-DATA 



NO 

1 

2 


1 

2 

3 

A 

5 

6 

7 

1 

A.AA23A 

9.52889 

#**«««** 

A.99669 

5.576A9 

6.18071 

6.80829 

7.A5810 

8.12897 

8.81966 

2 

9.52B89 

15.69856 


10.25533 

10,99757 

11.75A17 

12.52364 

13.30AAA 

1A.09A95 

14.89355 

3 

15.69856 

22.15508 


16.50825 

17.32085 

18.13A56 

18.94760 

19. 75806 

20.56406 

21.36371 

A 

22.15508 

27.95369 

««*«««** 

22.9362A 

23.70522 

2A.A6007 

25.19682 

25.91952 

26.62020 

27.29890 
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Illllll I ill 


5 

27.95369 

32.10503 


28.58264 

29.18385 

29.75542 

30.29550 

30.80228 

31.27395 

31.70877 

6 

32-10503 

33-59923 


32.46107 

32-77528 

33.04610 

33-27204 

33.45165 

33.58358 

33.66652 

7 

33.69923 

32-03130 


33.68056 

33.60943 

33.48485 

33-30590 

33.07176 

32-78170 

32.43505 

8 

32.03130 

26.71159 


31.56997 

31.05071 

30-47327 

29.83749 

29.14333 

28.39083 

27.50015 

9 

26.71159 

23-95033 

4 :^ 4 ^***** 

23.95033 







10 

23-95033 

21.18908 

* 4444:4441 

21- 18908 







11 

21- 18908 

17.94000 

4 ** 44*44 

20.61498 

20-08873 

19.61058 

19.18061 

18, 79879 

18.46497 

18.17885 

12 

17,94000 

18.45884 

*** 4 **** 

17.77343 

17.67878 

17.65372 

17.69541 

17. B0060 

17.96563 

18.18649 

13 

18.45884 

24.63803 

4 ** 44 *** 

19. 15903 

20.04390 

21.02410 

22.01216 

2^.92841 

23.70564 

24.29097 
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